Bulletin of Electrical Engineering and Informatics 
Vol. 11, No. 1, February 2022, pp. 9~19 
ISSN: 2302-9285, DOI: 10.1159 1/eei.v1 111.3074 o 9 


Two-port network model of wind turbine generator for three- 
phase unbalanced distribution system load flow analysis 


Rudy Gianto, Purwoharjono 


Department of Electrical Engineering, Tanjungpura University, Pontianak, Indonesia 


Article Info 


ABSTRACT 


Article history: 


Received May 19, 2021 
Revised Oct 12, 2021 
Accepted Jan 6, 2022 


Keywords: 


Induction generator 

Load flow 

Three-phase distribution system 
Two-port network 

Wind turbine generator 


An unbalanced electric power system can occur due to unequal lines and/or 
loads in each phase of the system. For a system with a high level of 
unbalance, the assumption that the system is balanced becomes incorrect. 
For a case like this, the system must be analyzed using three-phase 
approaches because single-phase analysis techniques can no longer be used 
to obtain accurate results. Therefore, the development of three-phase models 
of all components of the electric power system, including the wind turbine 
generating system (WTGS), is very important so that the system can 
properly be analyzed and its performance can correctly be evaluated. This 
paper proposes a method to incorporate WTGS in a three-phase load flow 
analysis of an unbalanced distribution system. The proposed method is based 
on the previously published single-phase two-port network model. In this 
paper, the single-phase two-port network model is modified and extended so 
that it can be applied to unbalanced systems. The model is then included in 


the three-phase distribution system load flow (DSLF) analysis. The proposed 
method is validated using 12-node and 19-node distribution systems. 
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1. INTRODUCTION 

It has been well acknowledged that penetration of wind turbine generating system (WTGS) in the 
electric power system will complicate the assessment and evaluation of the power system steady state and 
dynamic performances. This complication can be more profound in unbalanced electric power distribution 
systems since three-phase techniques have to be used in the analyses. An unbalanced electric power system 
generally occurs due to unequal lines and/or loads in each phase of the power system. For a system with a 
high level of unbalance, the assumption that the system is balanced becomes incorrect. For a case like this, 
the electric power system must be analyzed using three-phase approaches because single-phase analysis 
techniques can no longer be used to obtain accurate results [1]. 

To assess and evaluate the system steady state and dynamic performances of a power system 
embedded with WTGS, all of the system components (including the WTGS) need to properly be modeled 
[2]. In recent years, many investigations have been conducted in the context of modeling and integrating 
WTGS into steady state load flow analysis of an electric power system. For example, references [3]-[21] 
propose interesting methods to incorporate WTGS into load flow analysis. In [3]-[5], the induction generator 
equivalent circuit of the WTGS is converted to a three-node model. The three-node model application will 
modify the load flow analysis due to the addition of two branches and two PQ nodes to the original power 
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system where the WTGS is connected. With this method, a traditional load flow program can be used to 
solve the load flow problem, and modifications to the program source code are therefore not required. 

In [6]-[16], several steady-state mathematical models of fixed speed WTGS for load flow analyses 
have been developed. In the analyses, the equations of the WTGS models, combined with the existing power 
system equations (i.e., power system equations without the WTGS), are solved using some non-direct or 
iterative techniques. References [6]-[13] propose models which are formed based on the D or TI equivalent 
circuit of the WTGS induction generator. Reference [14] also proposes a steady-state model of fixed speed 
WTGS for load flow analysis. The model in [14] is derived based on the theory of a two-terminal network (or 
two-port network). In [15], [16] propose WTGS load flow models, which are obtained based on the Y or T 
equivalent circuit of the WTGS induction generator. The model in [15], [16] is derived using Kirchhoff and 
Ohm laws. In [17]-[21], DFIG-based variable speed WTGS models for steady state load flow analysis have 
been proposed. Equivalent circuits of the doubly fed induction generator (DFIG) have also been used in 
developing the models in [17]-[21]. 

All of the WTGS models proposed in [3]-[21] assume that the system was balanced so that single- 
phase load flow techniques can be used to analyze it. However, as mentioned earlier, the electric power 
distribution system is generally unbalanced and must be analyzed using three-phase load flow techniques. 
Against the above background, this paper proposes a three-phase WTGS model that can be applied to 
unbalanced distribution systems. The proposed method is based on the two-port network model, which has 
successfully been applied to a single-phase system [14]. Therefore, the main contributions of the present 
work are as shown: i) development of WTGS steady state model to be used in three-phase load flow analysis, 
and ii) modification and extension of the single-phase two-port network model so that it can facilitate three- 
phase load flow analysis under unbalanced system conditions. Moreover, validation studies using 12-node 
and 19-node distribution systems are also presented and discussed in this paper. To be systematic, the paper 
is organized as: section 2 discusses the formulation of the three-phase distribution system load flow (DSLF) 
problem. Section 3 addresses the single-phase two-port network model of WTGS. Extension of the single- 
phase two-port network model to a three-phase two-port network model is given in section 4. Section 5 
discusses a study case where validation of the proposed method is presented and investigated. Finally, some 
essential conclusions of the paper are pointed out in section 6. 


2. FORMULATION OF THREE-PHASE DSLF PROBLEM 
It can be shown that the formulation of three-phase load flow problem of unbalanced distribution 
system will have the form in (1): 


{[diag (V%®°)] 1 (sare — sabe) me yabcyabc =0 (1) 
where V”* is vector of nodal voltage, Sg“ is vector of power entering the node, Sz“° is vector of power 


leaving the node, and ¥** is system admittance matrix. 
For n-node distribution system, the formulations for V°, S°, S1% and Y*”¢ are given by: 


a. Save Spe yape yebe p yee 
yore = ngabe = Sue = yobe = Ya yee ; yee (2) 
yabe sge save vane ype : yabe 


where: 


b 
vi Sh yee Ya Ya WP YEP YS 
a sap = si ijr- s&h = yg væ Yge =|vie ve W| © 
Shi yea yg? yc ye sei yee 


J y 


It can be seen that (1) is a set of nonlinear equations which has to be solved in the load flow analysis 
to evaluate the steady state condition of the distribution system. As it is nonlinear, iterative techniques, for 
examples: Newton-Raphson or trust-region methods [22], [23], can be used to find the solution to the set of 
equations. 
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3. TWO-PORT NETWORK MODEL OF WTGS 
3.1. WTGS structure and two-port network 

A fixed-speed WTGS connected to bus k of an electric power distribution system is shown in 
Figure 1. In fixed-speed WTGS, a squirrel cage induction generator (SCIG) is usually used to convert the 
turbine mechanical power to electrical power. In Figure 1, Sz is WTGS electrical power output and Pm is 
turbine mechanical power input. Steady state equivalent circuits of the SCIG is given in Figures 2(a) and 2(b) 
[3]-[6]. In Figure 2(a), Rs and Xs are stator resistance and reactance, Rr and Xp are are rotor resistance and 
reactance, Re and Xm are core resistance and reactance, Vs and Js are stator voltage and current, Vr and Jp are 
rotor voltage and current, and s is induction generator slip. Also the impedance values of Zs, Zr and Zm in 
Figure 2(b) can be calculated using: 


Zs = Rs + jXs (4a) 
Zr = Ret jXp (4b) 
Zu = JR-Xm/(Re + jXm) (4c) 


Electric 
Power 
Distribution 
System 


kee DE RrtjXr 
@ 
Re JXm Rp(1-s)/s = 
& 
(a) 


Figure 2. Steady state equivalent circuit of SCIG (a) in terms of resistances and reactances and (b) in terms of 
impedances 


Figure 3 shows a general two-port network containing passive impedances/admittances [14]. In 
Figure 3, A, B, C and D are the general two-port network constants. The values of these constants depend on 
the impedances/admittances of the network components. For the two-port network in Figure 3, it can be 
shown that the voltage/current relationships are of the forms: 


1 


V, — EOM 
R — AD-BC 


(DV; + BI;) (5a) 
1 


Ip = —— 
R — AD-BC 


(CV; + Als) (5b) 

The two-port network model of WTGS is obtained by viewing the SCIG equivalent circuit in 
Figure 2(b) as a two-port network. For the equivalent circuit in Figure 2(b), the formulas for the 4, B, C and 
D constants of the two-port network will have the following forms [14]: 


A=14+2Zs/Zy (6a) 
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B =Z; + Zr +ZsZr/Zu (6b) 
C=1/Zy (6c) 
D=1+4+2Zp/Zy (6d) 


where Zs, Zr and Zm are calculated using (4). Development of the WTGS model using the two-port network 
theory will be discussed in the next section. 


Figure 3. General two-port network 


3.2. Two-port network model of WTGS 

Derivation of the two-port network model of WTGS is based on WTGS power formulas, i.e. WTGS 
electrical power output and turbine mechanical power input. By looking at Figure 2, the formulas for WTGS 
electrical power output (S,) and turbine mechanical power input (Pm) can be written as: 


Sg =P, + jQ = VsI (7a) 
and: 
Fn = VriR (7b) 


Substituting (5) into (7b) gives: 


Pa = |r Os + Bls)] [eae (CW + Ao] (8) 


On using (7a) in (8) and rearranging, can be obtained: 


Pn = ZV, + Z2Sg + Z3S + ZuSy (9) 
where: 
Dc* DA* * 
45 GD-RDjaD rO A ~ (AD-BC)(AD-BC)*’ Va = Vols (10a) 
= BC* on BA* niin = SRS 
23 = (acaBe)(ab—pey lA ~ (AD-BC)(AD-BC)*’"" — V,ve (10b) 


As shown in (9) is the steady state model of WTGS for single-phase load flow analysis. Extension 
of the model to three-phase load flow analysis will be discussed in the next section. 


4. THREE-PHASE TWO-PORT NETWORK MODEL OF WTGS 
For three-phase system, in (9) can be modified to become: 


pabe = Zabey abe + Sar + Zare( sare). + 700e gabe (11) 
where: 
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s$ ACN RA ze o 
Sart = S? ygce = VP (V) ; gabe en 9 ( 
Sg VP VPY 


Pn/3 
pope = [es 
Pn/3 


Diagonal elements of Z/*° matrix in (12) will have the following forms: 


Z? E DPCP* . 7P? = DP AP* 

1 (APDP-BPCP)(APDP-BPCPC)*?2 (APDP—BPCP)(APDP—BPCPC)* 
zP BPcP* . zP _ BP AP* 

3 (APDP-BPCP)(APDP-BPCPC)*’" "4  (APDP-BPCP)(APDP-BPCPC)* 


where: p=a, b, c; and Æ, BP, C’ dan D? are calculated using: 
AP=14+Z0/Z4, 
B? =Z? + ZP + Z0Z0/Zyy 
C? = 1/77, 
D? =1+Zh/Zy 

Based on (4), the formulations for Zs’, Zp’, and Zy” in (14) are: 
Ze = RE + jX? 
Ze = RR +jXR 


Zi = JREXE /(RE + JX) 


(13a) 


(13b) 


(14a) 
(14b) 
(14c) 


(14d) 


(15a) 
(15b) 


(15c) 


13 


To carry out three-phase load flow analysis for distribution system embedded with WTGS, the first 


step is to rearrange (11) to the following form: 


abc abcyabc abc gabe abc (cabc )* abc çabc _ 
pabe _ gabeyabe _ ggbc gabe _ Zebe( gabe)" _ abe gabe = 0 


(16) 


Then, for the system nodes containing WTGS, (16) is combined with (1). These combined 
equations, together with the remaining nodal equations, are solved simultaneously. Table 1 shows all of the 
quantities and equations in the complete formulation of a three-phase load flow problem for an electric power 


distribution system containing WTGS. 


It should be noted that in (16), the voltage at WTGS node or WTGS terminal is also the voltage of 


SCIG stator, and the power injection at WTGS node is also the power output of WTGS, therefore: 


abc _ yabc 
vere = ye 
abc — çabc 
SGj = S3 


where j is the node where the WTGS is connected. 


(17a) 


(17b) 


Table 1. Known and unknown variables in three-phase DSLF formulation for system with WTGS 


Node Equation Known variable Unknown variable 
Substation a) yare, gabe. ye sare 
Load (1) ye, sabe. gabe yo 
WIGS = (1) & (16) ¥9°; SBS; PAP; ZE; ZR; Zh Vere Sane 
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5. CASE STUDY 
5.1. Test system 

In this section, application of the proposed method in DSLF analysis is investigated. The test 
systems that are used in the investigation are described as shown in: 

- System 1 is an 11 kV distribution network with the total three-phase load of 4,305 kW and 1,215 kVAR 
[22]. Single line diagram of this 12-node balanced distribution system is shown in Figure 4. Since the 
system is balanced, it can actually be analysed using single-phase approach. However, to validate the 
proposed method, three-phase approach will be used instead. Tables 2 and 3 show the data of system 1. 
WTGS with power rating of 3,000 kW is assumed to be connected to node 12. 


Figure 4. System 1 


Table 2. Data of system 1 


Line No. Sending bus _ Receiving bus R* X* Pist O,** 
1 1 2 1.093 0.455 60 60 
2 2 3 1.184 0.494 140 30 
3 3 4 2.095 0.873 155 55 
4 4 5 3.188 1.329 130 30 
5 5 6 1.093 0.455 120 15 
6 6 T 1.002 0.417 155 55 
7 7 8 4.403 1.215 145 45 
8 8 9 5.642 1.597 140 40 
9 9 10 2.890 0.818 135 30 
10 10 11 1.514 0.428 140 30 
11 11 12 1.238 0.351 115 15 


*Line resistance/reactance in ohm/phase 
**Load connected to receiving bus in kW and kVAR (perphase) 


Table 3. WTGS induction generator data for system 1 


Parameter Resistance/reactance value (Ohm) 
Stator circuit R RŻ=R=1.21 
XT X}=X=6.05 
Rotor circuit RORS=R,=1.21 
XX =X,=6.05 
Magnetic core circuit RORP=RS=12,100 
Xn Xo =Xn=605 


- System 2 is a 19-node unbalanced distribution system adopted from [24], [25]. Single line diagram of 
the system is shown in Figure 5. It is an 11 kV distribution network with the total load of 442.75 kW 
and 152.45 kVAR (Phase-a: 151.06 kW and 52.66 kVAR; Phase-b: 143.07 kW and 48.45 kVAR; 
Phase-c: 148.62 kW and 51.34 KVAR. The system line and load data can be found in Tables 4 dan 5. 
WTGS with power rating of 300 kW is assumed to be connected to node 19. Data of the WTGS 
induction generator is given in Table 6. 


Figure 5. System 2 
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Table 4. Test system line data 


Z 
ô 
m 
5 
© 


Self impedance (Ohm) 


Mutual impedance (Ohm) 


SCIDAKRWNHH 
T 
EN 


3.0x(1.56090+j0.67155) 
5.0x(1.56090+j0.67155) 
1.5x(1.56090+j0.67155) 
1.5x(1.56090+j0.67155) 
1.0x(1.56090+j0.67155) 
2.0x(1.56090+j0.67155) 
2.5x(1.56090+j0.67155) 
3.0x(1.56090+j0.67155) 


3.0x(0.52030+j0.22385) 
5.0x(0.52030+j0.22385) 
1.5x(0.52030+j0.22385) 
1.5x(0.52030+j0.22385) 
1.0x(0.52030+j0.22385) 
2.0x(0.52030+j0.22385) 
2.5x(0.52030+j0.22385) 
3.0x(0.52030+j0.22385) 


9 9-10 5.0x(1.56090+j0.67155) 5.0x(0.52030+j0.22385) 
10 10-11 1.5x(1.56090+j0.67155) 1.5x(0.52030+j0.22385) 
11 10-12 1.5x(1.56090+j0.67155) 1.5x(0.52030+j0.22385) 
12 11-13 5.0x(1.56090+j0.67155) 5.0x(0.52030+j0.22385) 
13 11-14 1.0x(1.56090+j0.67155) 1.0x(0.52030+j0.22385) 
14 12-15 5.0x(1.56090+j0.67155) 5.0x(0.52030+j0.22385) 
15 12-16  6.0x(1.56090+j0.67155) 6.0x(0.52030+j0.22385) 
16 14-17 3.5x(1.56090+j0.67155) 3.5x(0.52030+j0.22385) 
17 14-18  4.0x(1.56090+j0.67155) 4.0x(0.52030+j0.22385) 
18 15-19 4.0x(1.56090+j0.67155) 4.0x(0.52030+j0.22385) 
Table 5. Test system load data 
Node Phase-a Phase-b Phase-c 
P(kW) Q(kVAR) P(kW) Q(KVAR) P(kW) Q(KVAR) 
1 0 0 0 0 0 0 
2 10.38 5.01 5.19 2:52 10.38 5.01 
3 11.01 5.34 5.19 252 9.72 4.71 
4 4.05 1.95 5.67 2.76 6.48 3.15 
5 6.48 3.15 5.19 2:52. 4.53 2.19 
6 4.20 2.04 3.09 1.50 2.91 1.41 
7 9.72 4.71 8.10 3.93 8.10 3.93 
8 7.44 3.60 5.34 2.58 3.39 1.65 
9 12.30 5.97 14.91 7.23 13.29 6.42 
0 3.39 1.65 4.20 2.04 2.58 1.26 
1 7.44 3.60 7.44 3.60 1.01 5.34 
2 9.72 4.71 8.10 3.93 8.10 3.93 
3 4.38 2.13 5.34 2.58 6.48 3.15 
4 3.09 1.50 3.09 1.50 4.05 1.95 
5 14.38 2.13 14.86 2.34 6.96 3.36 
6 17.77 3.78 20.38 5.01 7.77 3.78 
7 16.48 3.15 14.86 2.34 4.86 2.34 
8 15.34 2.58 15.34 2.58 5.52 2.67 
9 18.76 4.23 20.05 4.86 7.14 3.45 


Table 6. WTGS induction generator data for system 2 


Parameter Resistance/reactance value (Ohm) 
R R}=R=0.080667 


Stator circuit 


Rotor circuit 


Magnetic core circuit 


XTX =X=0.403333 

R“RL=R,=0.080667 

XXP=X,=0.403333 
R&RL=R-=20166.6665 
Xn Xn? =Xn=1008.3335 


5.2. Results and discussion (system 1) 

Results of the load flow calculation of system 1 are given in Tables 7 and 8. WTGS model, as 
proposed in section 4, is used in the calculation. It is to be noted that since the system is balanced, a single- 
phase approach such as the multiple-node model proposed in [3]-[5] can also be applied to analyze the 
system. For the purpose of comparison, results obtained by the multiple-node model are also given in the 
present work and are shown in Table 9. It can be seen that the results of the proposed method (Table 7 and 
Table 8) and those of the multiple-node method (Table 9) are in exact agreement, which indicates the validity 


of the method proposed. 
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Table 7. WTGS terminal voltage magnitude (system 1) 
Pn (KW) Phase-a (kV) Phase-b (kV) Phase-c (kV) 


300 9.86 9.86 9.86 
600 10.13 10.13 10.13 
900 10.38 10.38 10.38 
1200 10.62 10.62 10.62 
1500 10.84 10.84 10.84 
1800 11.05 11.05 11.05 
2100 11.25 11.25 11.25 
2400 11.44 11.44 11.44 
2700 11.62 11.62 11.62 
3000 11.80 11.80 11.80 


Table 8. Output power of WTGS (system 1) 
Pn (kW) Phase-a (kW, kVAR) Phase-b (kW, kVAR) Phase-c (kW, kVAR) 


300 91.55-j160.67 91.55-j160.67 91.55-j160.67 
600 190.39-j173.26 190.39-j173.26 190.39-j173.26 
900 288.87-j 187.36 288.87-j187.36 288.87-j187.36 
1200 387.03-j202.79 387.03-j202.79 387.03-j202.79 
1500 484.91-j219.38 484.91-j219.38 484.91-j219.38 
1800 582.53-j237.05 582.53-j237.05 582.53-j237.05 
2100 679.93-j255.68 679.93-j255.68 679.93-j255.68 
2400 777.11-j275.20 777.11-j275.20 777.11-j275.20 
2700 874.10-j295.56 874.10-j295.56 874.10-j295.56 
3000 970.90-j3 16.69 970.90-j3 16.69 970.90-j3 16.69 


Table 9. Results obtained by multiple-node model [3]-[5] 
Pa (kW) WTGS voltage (kV) WTGS output (kW, KVAR) 


300 9.86 91.55-j160.67 
600 10.13 190.39-j173.26 
900 10.38 288.87-j 187.36 
1200 10.62 387.03-j202.79 
1500 10.84 484.91-j219.38 
1800 11.05 582.53-j237.05 
2100 11.25 679.93-j255.68 
2400 11.44 777.11-j275.20 
2700 11.62 874.10-j295.56 
3000 11.80 970.90-j316.69 


5.3. Results and discussion (system 2) 

Tables 10-13 show the results of load flow analysis for the 19-node distribution system in Figure 5. 
To make a better observation, the results are also presented in graphical forms (see Figures 6-9). Due to the 
system unbalance, it can be seen from the tables that the electrical quantities in the three phases are not the 
same. Table 10 shows that with the increase in turbine mechanical power and WTGS output, the voltage 
profile improves (see also Figure 6). On the other hand, results in Table 11 show that as the turbine 
mechanical power input increases, the active power output of the WTGS also increases linearly (see also 
Figure 7). However, the WTGS power output is slightly smaller than the turbine mechanical power due to 
power loss in the WTGS. Table 11 also shows that the reactive power output of WTGS is negative, which 
means that the WTGS absorbs reactive power from the system. This reactive power is used by the WTGS 
induction generator for the purposes of magnetizing the core circuit. Figure 8 shows the variation of WTGS 
reactive-power demand against turbine mechanical power. It can be seen that the relationship is also linear, 
which indicates that the demand for reactive power increases as the WTGS active power output increases. 


Table 10. WTGS terminal voltage magnitude (system 2) 
Pan (kW) Phase-a (kV) Phase-b (kV) Phase-c (kV) 


30 10.0966 10.0749 10.0838 
60 10.1797 10.1587 10.1674 
90 10.2615 10.2410 10.2494 
120 10.3417 10.3219 10.3299 
150 10.4205 10.4012 10.4090 
180 10.4981 10.4792 10.4869 
210 10.5744 10.5559 10.5633 
240 10.6494 10.6315 10.6387 
270 10.7234 10.7059 10.7128 
300 10.7962 10.7790 10.7858 
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Table 11. Output power of WTGS (system 2) 
Pm(kW)  Phase-a (kW, kVAR) Phase-b (kW, KVAR) Phase-c (kW, KVAR)  3-Phase (kW, KVAR) 


30 8.31-j33.69 8.32-j33.55 8.32-j33.60 24.95-j100.84 
60 18.28-j34.25 18.29-j34.11 18.29-j34.17 54.86-j102.53 
90 28.25-j34.82 28.26-j34.68 28.26-j34.74 84.77-j104.24 
120 38.22-j35.38 38.23-j35.25 38.23-j35.30 114.68-j105.93 
150 48.19-j35.94 48.20-j35.81 48.20-j35.86 144.59-j107.61 
180 58.16-j36.50 58.17-j36.37 58.16-j36.43 174.49-j109.30 
210 68.13-j37.06 68.14-j36.94 68.13-j36.99 204.40-j110.99 
240 78.10-j37.62 78.10-j37.50 78.10-j37.55 234.30-j112.67 
270 88.06-j38.18 88.07-j38.06 88.07-j38.10 264.20-j114.34 
300 98.03-j38.74 98.04-j38.61 98.03-j38.66 294.10-j116.01 


Table 12. Substation power (system 2) 
Pn(kW)  Phase-a(kW, KVAR) Phase-b (kW, KVAR) Phase-c (KW, KVAR) 3-Phase (kW, KVAR) 


30 179.69+j99.95 169.824494.93 176.84+j98.43 526.35+j293.32 
60 168.40+j99.95 158.49+j94.92 165.52+j98.41 492.41+j293.28 
90 157.30+j 100.03 147.354j94.99 154.40+j98.48 459.05+j293.50 
120 146.38+j100.19 136.40+j95.14 143.46+j98.62 426.24+j293.95 
150 135.64+j 100.42 125.624j95.36 132.70+j98.84 393.96+j294.61 
180 125.06+j 100.72 115.02+j95.66 122.10+j99.12 362.19+j295.50 
210 114.65+j101.08 104.58+j96.02 111.67+j99.48 330.90+j296.58 
240 104.39+j101.52 94.29+)96.45 101.40+j99.90 300.08+j297.87 
270 94.28+j102.01 84.15+j96.94 91.27+j100.39 269.70+j299.34 
300 84.31+j102.57 74.16+j97.49 81.28+j100.93 239.75+j300.99 


Table 13. Distribution line loss (system 2) 
Pm (kW) Phase-a (kW, kVAR) Phase-b (kW, kVAR) Phase-c (kW, KVAR) 3-Phase (kW, kVAR) 


30 36.944) 13.60 35.07+j12.93 36.54+j13.49 108.55+j40.02 
60 35.62+j13.04 33.71+j12.36 35.19+j12.90 104.52+j38.30 
90 34.49+j12.55 32.54+j11.86 34.04+j12.40 101.07+j36.81 
120 33.54+j12.15 31.56+j11.44 33.07+j11.98 98.17+j35.57 
150 32.77+j11.82 30.75+j11.10 32.28+j11.64 95.80+j34.56 
180 32.16+j11.56 30.12+j10.84 31.64+j11.35 93.92+j33.75 
210 31.72+j11.36 29.65+j10.63 31.18+j11.15 92.55+j33.14 
240 31.43+j11.24 29.32+j10.50 30.88+j11.01 91.63+j32.75 
270 31.28+j11.17 29.15+j10.43 30.72+j10.95 91.15+j32.55 
300 31.28+j11.17 29.13+j10.43 30.69+j10.93 91.10+j32.53 
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Figure 6. Variation of voltage magnitude Figure 7. Variation of WTGS active power 


Table 12 shows that with the WTGS installed in the system, besides the voltage profile 
improvement, another advantage is that the power supply from the distribution system substation can be 
reduced since some of the loads can be fed by the WTGS (see also Figure 9). Power supply from distribution 
system substation usually comes from the conventional generator with a non-renewable energy source (fossil 
fuel) that is not environmentally friendly. Moreover, as the WTGS active power output increases, the line 
power loss decreases (see Table 13). This line loss decrement is due to the voltage profile improvement as 
the WTGS active power output increases. The above results confirm the validity of the method proposed. 
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Further validation can be obtained by examining the results where the substation output plus WTGS output is 
always equal to the total system load plus total line loss. It is to be noted that the line loss has been calculated 
based on the impedances and currents of the distribution lines. 
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Figure 8. Variation of WTGS reactive power Figure 9. Variation of substation power 


6. CONCLUSION 

Most of the previous WTGS models assumed that the system was balanced, and single-phase load 
flow approaches were used to analyze the system. However, electric power distribution systems are generally 
unbalanced and must be analyzed using three-phase load flow techniques. Therefore, the development of a 
three-phase model of all components of the electric power system (including WTGS) is very important so 
that the system can be analyzed and its performance can properly be evaluated. A simple method to 
incorporate WTGS in load flow analysis of a three-phase unbalanced distribution system has been proposed 
in this paper. The proposed method is based on a previously published single-phase two-port network model. 
In this paper, the single-phase two-port network model has been modified and extended so that it can be 
applied to three-phase systems. The model is then included in a three-phase load flow analysis of an 
unbalanced distribution network. An extensive case study has also been carried out in this paper. In the case 
study, the application of the proposed method to 12-node and 19-node distribution systems has been 
investigated. The investigation results confirm and verify the validity of the method proposed. 
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